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Abstract 

We use inelastic neutron scattering to demonstrate that at low tempera- 
tures, the low frequency magnetic fluctuations in YBa2Cu306.6 {Tc = 62.7 K) 
are incommensurate, being found at positions displaced by ^5 (0.057 it 0.006 
r.l.u.) along the [vr, vr] direction from the wave vector (vr, vr) associated with 
the antiferromagnetic order of the parent insulator, YBa2Cu306. The dy- 
namical susceptibility }("{(\,uj) at the incommensurate positions increases on 
cooling below T^, accompanied by a suppression of magnetic fluctuations at 

the commensurate points. 

PACS numbers: 74.72.Bk, 61.12.Ex 
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Knowledge of the spin dynamical properties of the cuprates are crucial to the under- 
standing of high-temperature (T^) superconductivity. An important issue is the symmetry 
of the imaginary part of the dynamical susceptibility, x"(q, cj), probed directly by neu- 
tron scattering. Over the past several years, intensive experimental work on the single 
layer La2_^Sr^Cu04 (214) and the bilayer YBasCusOr.^ [(123)07_^] cuprates 
has yielded valuable information concerning the magnetic response in the normal and super- 
conducting states. For the 214 family, magnetic fluctuations were found at incommensurate 
positions from the antiferromagnetic ordering position (7r,7r) [|l|]. For (123)07_2:! the situa- 
tion is more subtle. While Rossat-Mignod and coworkers only detected magnetic fluctuations 
at the commensurate position (vr, vr) Tranquada et al. noticed that the wave vector 
dependence of x"(q, tu) in a (123)06.6 (^c = 53 K) crystal was better described by a pair 
of identical Gaussians displaced symmetrically from the (vr, vr) point, suggesting possible in- 
commensurate fluctuations. In a subsequent experiment, Sternlieb et al. 0] have suggested 
that the wave vector dependence of the susceptibility is independent of energy for all ener- 
gies between 2 mV and 40 meV. Although the nontrivial wave vector dependence behavior 
was again observed, x"(q, u;) was not found to be different below and above T^. No flrm 
conclusions about the commensurability and symmetry of x"(q, cj) were reached 

In this Letter, we present inelastic neutron scattering data which resolves the issue of 
commensuration in (123)06.6- We show that the magnetic response is complex with incom- 
mensurate fluctuations for energies below the commensurate resonance ||^ at low tempera- 
tures. The low frequency spin fluctuations change from commensurate to incommensurate 
on cooling with the incommensuration flrst appearing at temperatures somewhat above T^. 
Lowering the temperature suppresses the spin fluctuations at the commensurate points, 
accompanied by an increase in the susceptibility at the incommensurate positions. 

Our results shed light on several long standing theoretical predictions. If the (123)07_a, 
system is indeed a d-wave superconductor, ci-wave gap nodes could yield incommensurate 
peaks below |]TD|]. In this scenario, the scattering of (123)06.6 is expected to be commensu- 



rate in the normal state and incommensurate in the superconducting state |Tl]]. No neutron 



scattering experiment, until now, has observed this behavior. Second, even though the differ- 
ences in the commensurabihty in (123)07_i. and 214 can be explained from the Fermi surface 
differences between the two compounds [O, such models have not been entirely reconciled 



with the nuclear magnetic resonance (NMR) measurements The reconciliation of 

the neutron and NMR results requires detailed information about the structure of x"(q, cu) 
||T5| . The third theoretical conjecture stems from recent neutron scattering experiments of 
Tranquada et al. [|l^] which suggest that the incommensuration in superconducting 214 may 
be associated with a spatial segregation of charge or charge density wave correlations |]17 



If the idea of dynamical microphase separation in the Cu02 plane asserted by Emery and 
Kivelson |1^ is relevant for the high-T^ superconductivity, one would expect incommensu- 



rate spin ffuctuations in other cuprate superconductors like Bi2Sr2CaCu208-x (2212) and 
(123)07_a;. Indeed, recent neutron scattering experiments by Mook and Chakoumakos p!9[ 



show that the spin ffuctuations in 2212 are also incommensurate. Thus, incommensurability 
may be a common feature for all cuprate superconductors. 

The neutron scattering measurements were made at the High-Flux Isotope Reactor at 
Oak Ridge National Laboratory using the HB-1 and HB-3 triple-axis spectrometers. The 
fabrication and characteristics of our single- crystal sample of (123)06.6 (weight 25.59 grams 
and Tc = 62.7 K) were described in detail previously f^. The major difficulty in study- 
ing spin ffuctuations in the (123)07_2,. system is to separate the magnetic scattering from 
(single- and multi-) phonon and other spurious processes. While multi-phonon scattering 
usually has a simple wave vector dependence and spurious events such as accidental Bragg 
scattering can be identified by checking the desired inelastic scan in the two- axis mode [Q, 
two approaches can be used to separate magnetic from single-phonon scattering. The first 
approach is to perform neutron polarization analysis I^O] which, in principle, allows an un- 



ambiguous separation of magnetic and nuclear scattering. This method has been successfully 
employed to identify the magnetic origin of resonance peaks for ideally and underdoped 
(123)07_^. However, this advantage comes at a considerable cost in intensity which 
makes the technique impractical for observing small magnetic signals. The second approach 



is to utilize the differences in the temperature and wave vector dependence of the phonon 
and magnetic scattering cross sections. While phonon scattering gains intensity on warming 
due to the thermal population factor, the magnetic signal usually becomes weaker because it 
spreads throughout the energy and momentum space at high temperatures. Thus, in an un- 
polarized neutron measurement the net intensity gain above the multi-phonon background 
on cooling at appropriate wave vectors is likely to be magnetic in origin. 

Figure 1(a) depicts the reciprocal space probed in the experiment with a* (=l/a), b* 
(=1/6) directions shown in the square lattice notation. The momentum transfers {qx, Qy, Qz) 
in units of A"-*^ are at positions {H, K, L) = {qxa/2TT, qyb/27T, qzc/27r) reciprocal lattice units 
(r.l.u.). We first describe measurements made in the {H,H,L) zone. Our search for the 
magnetic fluctuations was done with the filter integration technique first developed to 



study the chain fluctuations in (123)06.93. This technique is excellent for isolating scatter- 
ing from lower dimensional objects and relies on integrating the energy along wave vector 
direction [0, 0, L] perpendicular to the scan direction [H, H, 0]. To estimate the energy inte- 
gration range of the technique, we note that the scattered intensity for acoustic modulations 
in (123)07_a; is proportional to the in-plane susceptibility x"{lx,%,^) 



k 1 
I{q,uj) oc -^|/cu(q)psin2(-A2;g^)[n(^) + l]x"(gx,gy,w), 

where fcj and kf are the initial and final neutron wave numbers, /cu(q) is the Cu^+ mag- 
netic form factor, Az (=3.342 A) the separation of the Cu02 bilayers, q the total momentum 
transfer (|qp = ql + qy + q"^), and [n{u) + 1] the Bose population factor. The sohd line in Fig. 
1(b) shows the calculated /(q, u) at (vr, vr) as a function of energy transfer (along g^) assum- 
ing x"{Qx,'ly,^^) = F{.(lxAy)x" il^) ^F{.(lxAy) 1^3 • Although there are two broad peaks 
in the figure, the observed intensity will mostly stem from fluctuations around the lower 
energy one (10 < AE < 30 meV) because of the decreased resolution volume at large energy 
transfers. Since room temperature triple-axis measurements show no detectable magnetic 
peaks at (tt, tt) below 40 meV (see Figs. 2 and 3), we have used the integrated scan at 295 
K as the background and assumed subsequent net intensity gains above the multi-phonon 



background at lower temperatures are magnetic in origin. Figure 1(c) shows the result at 
different temperatures. At 200 K, the magnetic ffuctuations are broadly peaked at (7r,7r). 
On cooling to 150 K and 100 K, the peak narrows in width and grows in intensity but is 
still well described by a single Gaussian centered at (tt, tt). At 65 K, the data show a ffattish 
top similar to previous observations 0] . Although detailed analysis suggests that the profile 
is better described by a pair of peaks (Lorentzian or Lorentzian-squared line-shape) than a 
single Gaussian, the most drastic change in the profile comes in the low temperature super- 
conducting state. Rather than the expected single peak, two peaks at positions displaced 
by ±5 (0.057 ± 0.006 r.l.u.) from H = 0.5 are observed, accompanied by a drop in the spin 
fluctuations at the commensurate position. The observation of sharp incommensurate peaks 
with the filter integration technique suggests that the incommensuration must be weakly 
energy dependent in the integration range. 

Although the integration technique is excellent in finding small peaks from the scattering 
of lower dimensional objects, it is important to confirm the result with conventional triple- 
axis measurements and to determine the symmetry of the incommensuration. For this 
purpose, we have realigned the sample in the {H, 3H, L) zone. If the 15 K profile in Fig. 1(c) 
stems from an incommensurate structure with peaks at (0.5±5, 0.5±(5) [see Fig. 1(a)], scans 
along the [H, 3H] direction are expected to peak a.t H = 0.477 and 0.523 r.l.u. for 6 = 0.057. 
On the other hand, if the underlying symmetry is identical to that of 214 [rotated 45° from 
Fig. 1(a)], the incommensuration in a [H,3H] scan should occur aX H = 0.466 and 0.534 



r.l.u. Figure 2 summarizes the result at 24 meV [^. The scattering at room temperature 
shows no well defined broad peak around (vr, vr), but at 70 K a two peak structure emerges. 
On cooling below Tc, the spectrum rearranges itself with a suppression of fluctuations at 
commensurate point accompanied by an increase in intensity at incommensurate positions. 
The wave vectors of the peaks in the [H, 3H] scan are consistent with incommensuration 
at (0.5 ± 5,0.5 ±5). It may also be possible to interpret the data with other structures, 
however, we will assume the symmetry of Fig. 1(a) until more precise measurements are 
made. 



In previous work, superconductivity was found to induce a strong enhancement in the 
x"(q, cj) at (tt, tt) for ideally and underdoped (123)07_j, at the resonance po- 

sitions. Although the intensity gain of the resonance below Tc is accompanied by a sup- 
pression of fluctuations at frequencies above it for the underdoped compounds ll,^, no 
constant-energy scan data are available at energies above the resonance. In light of the 
present result at 24 meV for the (123)06.6 sample which has a resonance at 34 meV it is 
important to collect data at these frequencies. Thus, we undertook additional measurements 
with improved resolution (collimation of 50"-40"-40"-120" ) in the hope of resolving possible 
incommensuration at high energies. Figures 3(a) and (b) suggest that the fluctuations at 
the resonance energy are commensurate above and below Tc with no appreciable change in 
width. For an energy above the resonance (42 meV), the scan is featureless at room temper- 
ature but shows a well defined peak centered at (tt, tt) at 75 K. Although superconductivity 
suppresses the magnetic fluctuations [see inset of Fig. 3(d)], the wave vector dependence of 
the line-shape cannot be conclusively determined due to the poor instrumental resolution 
at this energy. Unfortunately, further reduction in resolution volume is impractical due to 
a concomitant drop in the scattering intensities. 

Since the earlier polarized work |^ has shown that for (123)06.6 the 34 meV resonance is 
the dominante feature of x"(^) ^) (^) ^) ^^e low temperature superconducting state, it 
is important to compare the newly observed incommensurate peaks to the intensity gain of 
the resonance. Figure 4 shows the difference spectra between 15 K and 75 K at frequencies 
below and above the resonance. In the energy and temperature range of interest (15 K to 
75 K), the phonon scattering changes negligibly and the Bose population factor [n{uj) + 1] 
modifies the scattered intensity at high temperatures by only 3% at 24 meV and less at 
higher energies. Therefore, the difference spectra in the figure can be simply regarded as 
changes in the dynamical susceptibility. Inspection of Figs. 4(a) and (b) reveals that the 
susceptibility at the incommensurate positions increases on cooling from the normal to the 
superconducting state, accompanied by a suppression of fluctuations at the commensurate 
point. Comparison of Fig. 4(c) to Figs. 4(a) and (b) indicates that the net gain in intensity at 
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the incommensurate positions below Tc is much less than that of the resonance. For an energy 
transfer of 42 meV, the intensity drop appears uniform throughout the measured profile, 
however, instrumental resolution may mask any possible incommensurate features. Figure 
4(e) plots a summary of the triple-axis measurement in superconducting state. Although 
there are only two constant-energy scans for frequencies below the resonance, these data 
nevertheless confirm the result of the integrated technique. 

In conclusion, we have found that the low frequency magnetic fluctuations in a (123)06.6 
sample are incommensurate at low temperatures with incommensuration first appearing at 
temperatures above Tc- The dynamic susceptibility at incommensurate positions increases 
on cooling below Tc, accompanied by a suppression of magnetic fluctuations at the commen- 
surate point. 
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FIGURES 

FIG. 1. (a) Diagram of reciprocal space probed in the experiment. The dashed arrow indicates 
the scan direction with the integrated technique while the solid arrow represents the triple-axis 
measurements, (b) Calculated scattered intensity /(q, w) as a function of energy transfer. The 
effective energy integration range is mostly from 10 to 30 meV. Note kf is parallel to Qz in this 
technique, (c) Integrated measurements in which the data at 295 K are subtracted from 200 K, 
150 K, 100 K, 65 K, and 15 K. The data are normalized to the same monitor count. The solid 
lines in the 100 K, 150 K, and 200 K data are fits to single Gaussians and linear backgrounds. The 
solid lines in the 15 and 65 K data are two Lorentzian-squared peaks on linear backgrounds which 
best fit the data. 

FIG. 2. Triple-axis scans along (iJ,3if, 1.7) at 24 meV for (a) 295 K, (b) 70 K, (c) 58 K, 
and (d) 50 K. Data at 295 K were collected with HB-1 while other scans were taken using HB-3. 
Since room temperature measurements at 42 meV [Fig. 3(c)] show no peak at (vr, tt), the weak 
structures in (a) are most likely due to phonon and/or spurious processes. The horizontal bar 
shows the resolution along the scan direction and the vertical resolution is 0.14 A~^. The positions 
of incommensuration at i7 0.48 and 0.53 r.l.u. are indicated by the arrows. Solid lines in (b)-(d) 
are two Lorentzian-squared peaks on a linear background. The increased scattering at H > 0.6 
r.l.u. is due to phonons. 

FIG. 3. Constant-energy scans along {H,3H, 1.7) with energy transfer of 34 meV at (a) 75 K, 
and (b) 15 K. Identical scans at 42 meV at (c) 295 K (•), 75 K (o), and (d) 15 K (o). Inset (•) 
shows the temperature dependence of the scattering at (0.5,1.5,-1.7) for AE = 42 meV where the 
arrow indicates Tc. The multi-phonon background in the 295 K data has been scaled to the value 
at 75 K for clarity. Solid lines are Gaussian fits to the data. The q width of the resonance after 
deconvolving the resolution (horizontal bar) is 0.23 A~^. A similar result is obtained at 42 meV. 
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FIG. 4. Difference spectra along {H,3H,1.7) between low temperature (< Tc) and high tem- 
perature Tc + 12 K) at (a) 24 meV, (b) 27 meV, (c) 34 meV, and (d) 42 meV. All data were 
taken with the same monitor units. To a very good approximation, the data can be regarded as 
the difference in susceptibility between 15 K and 75 K, i.e., x"(15 K) — x"(75 K). Solid lines are 
guides to the eye. (e) Summary of triple-axis measurements. Open squares indicate incommensu- 
rate positions. Solid and open circles are the resonance and fluctuations at 42 meV, respectively. 
The error bars show the energy resolution and the intrinsic q width (FWHM). 
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